Two-channel Kondo tunneling in triple quantum dot 
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The effective spin Hamiltonian of a triple quantum dot with odd electron occupation weakly 
connected in series with left (/) and right (r) metal leads is composed of two-channel exchange 
and co-tunneling terms. Renormalization group equations for the corresponding three exchange 
constants J;, J r and J; r are solved (to third order). Since J; r is relevant, the system is mapped 
on an anisotropic two-channel Kondo problem. The structure of the conductance as function of 
temperature and gate voltage implies that in the weak and intermediate coupling regimes, two- 
channel Kondo physics persists at temperatures as low as several Tk- At even electron occupation, 
the number of channels equals twice the spin of the triple dot (hence it is a fully screened impurity) . 

PACS numbers: 72.10.-d, 72.15.-v, 73.63.-b 



Motivation: In the present work, a simple configura- 
tion of localized moment in nanostructures is studied, 
where the two-channel Kondo Hamiltonian appears in 
resonance tunneling. Concrete experiment is proposed in 
order to elucidate the pertinent physics at T > Tk (the 
Kondo temperature). In the strong coupling regime, a 
multichannel Kondo system is known to be a non- Fermi 
liquid yj, but construction of simple theoretical mod- 
els pertaining to experimentally feasible setups is noto- 
riously elusive. Examples are magnetic impurity scat- 
tering, physics of two-level systems and Kondo lattices 
(see |2j, LJ for review) . Recent attempts to realize two- 
channel Kondo effect in tunneling through quantum dots 
using peculiar setups still await experimental manifes- 
tation. The problem is exemplified in tunneling through 
a simple quantum dot sandwiched between two metallic 
"left" (I) and "right" (r) leads. Starting from the single 
impurity Anderson model, it is tempting to think of the 
two leads as a source of two tunneling channels. However, 
if the two fermion lead operators c kaa {k = momentum, 
a = spin projection and a = l,r the lead index) are cou- 
pled to a single dot electron operator d a , one channel can 
always be eliminated by an appropriate rotation in I — r 
space 0: 

cucri = cos6 k c k(T i + sm0 fe Cfeo-2, 

Char = -Sm6 k C k crl + COs6 k C k<y2 - (1) 

As a result, only the standing wave fermions c ka \ con- 
tribute to tunneling. It is therefore natural to expect 
that a generic two-channel Hamiltonian can be realized 
only when the rotation Q cannot eliminate the second 
channel. Such situation may arise, e.g., when, instead of 
a simple quantum dot one has a nanoobject consisting of 
several dots so that c ka i and c kar are coupled to different 
dot electron operators. 

Model Hamiltonian: Consider a triple quantum dot 
(TQD) which consists of three dots l,f,r, connected in 
series to left and right leads (see Fig. 1). The figure 
defines also tunneling amplitudes V a and hopping ampli- 
tudes W a between the side dots and the central one. The 
source of Kondo screening is a Coulomb blockade in the 
central dot /, which is supposed to have a smaller radius 



(and hence a larger capacitive energy) than the two side 
dots, namely, Qj^> Qi, r (cf. [6(). The tunneling Hamil- 
tonian has the form H t = J2 a =i,r Uka VaC ktra d aa + H.c, 
and the rotation JJJ does not eliminate the second tun- 
neling channel. Note that direct tunneling through the 
TQD is suppressed due to electron level mismatch and 
Coulomb blockade, so that only co-tunneling mechanism 
contributes to the current. To quantify these elementary 
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FIG. 1: Triple quantum dot in series (upper panel) and 
single-electron energy levels of each individual dot e a = 
e a — V r ga =bare energy minus gate voltage (lower panel). 

statements consider the case of TQD with three electrons. 
The Hamiltonian of the isolated TQD reads 

H d = H £aU ' 7 

a=l,r,f \ cr / 

+ Y, (W a dt a dfa + H.c), (2) 

a=L,r 

where n aa = d\ a d aa , with J^aa n °°- = 3 - Tnis Hamil- 
tonian can be easily diagonalized in the space of three- 
electron states | A) of the TQD. The lowest ones are clas- 
sified as a ground state doublet low-lying doublet 
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excitation \d,2) and quartet excitation \q). The single 
electron levels e a = e a — V ga are tuned by gate voltages 
V ga such that [3 a = W a /A a < 1 (A Q = e a - e f ). To 
order (3 2 , the corresponding energy levels E\ are 



E dl =£f + El 



E r - 
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(3) 



The states |g?i) and \d,2) include components with two 
electrons on I or r dots (responsible for co-tunneling). 
Coupling with the leads admixes three-particle states |A) 
of the dot, with two and four particle states |A) (for 
brevity we include only the formers, addition of the lat- 
ters is straightforward). It is then useful to introduce 
diagonal and number changing dot Hubbard operators 
X AA = |A)(A| and X XA = |A> <A| respectively. The 
Hamiltonian of the whole system (leads, TQD and tun- 
neling) then reads 



H 



kaa 



\ AA kaa 



A V x 



AA 



H.c. 



(4) 



Here -Hf- < e ka < ^ 



are lead electron energies with 
nearly identical bandwidths Di « D r = D\. Within 
linear response, the Fermi energies are the same on both 
sides, 6fi = f-Fr = £f = 0. The tunneling amplitudes 
V AX = 14 (Al^o-a | A) depend explicitly on the lead index a 
and on the respective 2-3 particle quantum numbers A, A. 
Elimination of one of the two channels by the rotation 
is then impossible: adding I or r electron to a given state 
| A) results in different states |A). 

RG procedure and spin Hamiltonian: Following a 
two-stage RG procedure the low-energy Kondo tunneling 
through TQD may be exposed. First, the bandwidth is 
continuously reduced from D\ to Dq < D\ (see |(j and 
references therein), and the energy levels © are renor- 
malized by eliminating high-energy charge excitations 
0. If at the end of this procedure E dl (D ) < -D /2 
and Ed 2 (D ) — Ed^Do) exceeds Tk (defined below), 
charge fluctuations are quenched and the doublet ground 
state \di a= ^ t i) becomes a localized moment. Following 
the Schrieffer- Wolff transformation the spin Hamiltonian 
reads 

H s = JiS ■ Sf + J r S ■ s r + Ji r S ■ (si r + s r i). (5) 
The spin 1/2 operator S acts on \di a =i t [) whereas the 
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lead electrons spin operators are s Q = ^ 
The presence of the left-right spin operator s; r 
Sfefc' c k t j,i' T i-iv c k'vr is responsible for co-tunneling current. 
Moreover, J a = 8~/ 2 \V a \ 2 /3(e F - e a ) > (where 7 = 

Jl-Kff+ff)), While J W = -WlV r PlPr/3(e F - Ef) 



so that \Jir/J a \ is of order (3if3 r <C 1. The Hamil- 
tonian J5J) then encodes a two-channel Kondo physics, 
where the leads serve as two independent channels and 
T K = max{T Kl ,T Kr }, T Ka = D a er D "/^. 

In the second stage of the RG procedure, a poor- 
man scaling technique is used to renormalize the ex- 
change constants by further reducing the band-width 
Do — ► D. The pertinent fixed points are then identi- 
fied as D — > Tk &}■ Unlike the situation encountered in 
the single-channel Kondo effect, third order diagrams in 
addition to the usual single-loop ones should be included 
(see Fig.5 in Ref. Q and Fig. 9 in Ref. 0). Below we 
use D as an energy unit, hence J a , V a , W a , V ga , s a , T, 
D etc. now become dimensionless. With a = l,r and 
a = r,l the three RG equations for J; , J r , J\ r are, 

' /J " =-{Jl + Jl) + Ja(Jl + Jl + 2Jil 



d\nD 

dJlr 



-J lr (J, + J r ) + J lr (jf + J 2 r + 2jf r ). (6) 



d\nD 

On the symmetry plane Ji = J r = J, equations (j6j) re- 
duce to a couple of RG equations for Ji 2 = J ± J\ r 

dJi/dhxD = -J 2 + Ji{Jl + j\) (i = 1, 2), (7) 

subject to Ji(D = 1) = Jio- These are the well-known 
equations for the anisotropic two-channel Kondo effect 0. 
With 4> t = (Ji + J 2 - I)/ Ju Ci = fc - (fiio and Li(x) = 
x — ln(l + Ci/x) — 2 In a;, the solution of the system is 



u 



L, 



h\D (i = l,2). 



(8) 



The scaling trajectories in the sector (J; > J r > 0, J; r = 
0) and in the symmetry plane with < Ji r < J are shown 
in Fig. |21 Although the fixed point (1/2, 1/2, 0) remains 
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FIG. 2: Scaling trajectories for two-channel Kondo effect in 
TQD. 

inaccessible if Ji r 7^ 0, one may approach it close enough 
starting from an initial condition Ji r o Jio,J r o- Real- 
ization of this inequality is a generic property of TQD in 
series shown in Fig.l. 

Conductance: According to general perturbative ex- 
pression for the dot conductance 10], its zero-bias 
anomaly is encoded in the third order term, 

G Jl[Ji(T) + J r (T)}, (G = ^-). (9) 



G 0) 



Here the temperature T replaces the bandwidth D in 
the solution (jHJ). Let us present a qualitative discus- 
sion of the conductance G[J a (T)] (or in an experimen- 
talist friendly form, G(V ga ,T)) based on the flow dia- 
gram |2 (Strictly speaking, the RG method and hence 
the discussion below, is mostly reliable in the weak cou- 
pling regime T > Tk)- Varying T implies moving on a 
curve [Ji(T), J r (T), Ji r (T)] in three dimensional param- 
eter space (Fig. EJ, and the corresponding values of the 
exchange parameters determine the conductance accord- 
ing to equation @ . Note that if, initially, Jiq = J r o = Jo 
the point will remain on a curve [J(T), J(T), Ji r {T)} lo- 
cated on the symmetry plane. By varying V ga it is possi- 
ble to tune the initial condition (Jio, J r o) from the highly 
asymmetric case Jiq 3> J r o to the fully symmetric case 
Jia = Jro- For a fixed value of Ji r o the conductance 
shoots up (logarithmically) at a certain temperature T* 
which decreases toward Tk with | Jiq — J r o| and Ji r $. The 
closer T* is to Tk, the closer is the behavior of the con- 
ductance to that expected in a generic two-channel situ- 
ation. Thus, although the isotropic two-channel Kondo 
physics is unachievable in the strong coupling limit, its 
precursor might show up in the intermediate coupling 
regime. 

The conductance G(V ga , T) as function of T for several 
values of V ga and the same value of Ji r o is displayed by 
the family of curves in Fig0 For G displayed in curve a, 
T* /Tk ~ 3 and for T > T* it is very similar to what is ex- 
pected in an isotropic two-channel system. Alternatively, 
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FIG. 3: Conductance G in units of Go as a function of tem- 
perature (r = T/Tk), at various gate voltages. The lines 
correspond to: (a) the symmetric case Ji = J r (V g i = Vgr), 
(b-d) Ji » J r , with V g i - V gr = 0.03, 0.06 and 0.09. At 
t — > oo all lines converge to the bare conductance. 

holding T and changing gate voltages V ga enables an ex- 
perimentalist to virtually cross the symmetry plane. This 
is equivalent to moving vertically downward on Fig. [3] 
At high temperature the curves almost coalesce and the 
conductance is virtually flat. At low temperature (still 
above Tk) the conductance exhibits a sharp minimum. 
This is summarized in FigQ] 

Higher degeneracy and dynamical symmetries: 

The spin Hamiltonian J5J is expressible in terms of the 
dot (ground-state) spin S = 1/2 components Si (i = 
x, y, z), the generators of the group SU (2). Due to com- 
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FIG. 4: Conductance G in units of Go as a function of gate 
voltage at various temperatures (at the origin Ji — J r ). 

plete spin degeneracy Ed^ = E dl i, the dot Hamilto- 
nian itself is of course invariant under SU (2) but the 
hybridization leading to the exchange terms J a S • s Q 
clearly breaks it, allowing for dot spin-flips. In com- 
posite quantum dots, the degeneracy of spin states at 
the end of the first stage RG procedure might be much 
richer Q (and greater than 2). Our analysis shows that 
for the present model of TQD with N = 3 electrons 
there exists a scenario of level degeneracy in which the 
renormalized energies of the two doublets and the quar- 
tet are degenerate at the Schrieffer- Wolff limit, that is, 
E dl (D ) = E d2 (D ) = E q {D ) < -D /2. The corre- 
sponding wave functions |A) are vector sums of states 
composed of a "passive" electron sitting in the central 
dot and singlet /triplet (S/T) two-electron states in the 
I, r dots. Using certain combinations of dot Hubbard op- 
erators | A) (A' | one can now define two vector operators 
S and M such that S is the dot spin 1 operator respon- 
sible for transitions within the triplet states, while M 
accounts for S/T transitionsfgj- The operators S and M 
together with the spin 1/2 operator s of the central dot 
electron generate the group SO (A) x SU(2) specified by 
the Casimir operator S 2 + M 2 + s 2 — -y. The corre- 
sponding spin Hamiltonian, 

H * = J2 [ J a TS + J f TM ] ' S a + JtrS ■ (S lr + Br,), (10) 
a— l,r 

expressible in terms of its generators breaks that sym- 
m etry. Here Jj = 4 7 |1/ a | 2 /3(eF - e,) > 0, J a 5T = 

y/l-±(j3?+/%)JZ and J lr = -iVtVrftPr/iep - e f ), 

so that \Ji r \/J T w /3 2 . Both S and M vectors are in- 
volved in Kondo screening, a situation much richer than 
the one described by the single impurity Hamiltonian J^} • 
As far as the Kondo physics is concerned, the number of 
channels n = 2 exactly equals twice the spin, 2S. There- 
fore, the physics in the strong coupling limit is similar 
to that of the single channel one. The two stable fixed 
points are (Jj, jf T , Jj, jf T ) = (00,00, 0,0) implying 
T K = max{T Kl ,T Kr }, T Ka = {exp [-l/( jj + J a ST )] } . 
In the weak coupling limit, systems with 2S > n and 
2S = n might have different physics. 
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Finally, let us discuss the experimentally relevant sit- 
uation of changing electron occupation iV in the same 
TQD device. For N = 2, the lowest-energy states con- 
sist of two singlets \S a ) and two triplets \T a ) (one electron 
in the left and right dot and one electron in the central 
dot). The next charged transfer excitons consist of a sin- 
glet \Sir) and a triplet \T; r ) (one electron in the left and 
right dots). The corresponding energy levels are 

2W 2 

E Sa = e f +e a - a - WaPa, 
E T a = e f + e a - W a f3a, 

E Slr = E Tlr =ei+s r + Wi(3i + W r (3 r . (11) 

where a = W a /A a < 1, (A a = e„ - £/). 

The most symmetric case is realized when si = e r = e 
and Vi = V r = V. If at the end of the first stage RG 
procedure one has Es, ~ £r, ~ E,s r ~ Er r , the TQD 
possesses a P x 5*0(4) x SO (A) symmetry (with P is 
I <-> r permutation). Following an RG procedure and 
a Schrieffer- Wolff transformation, the spin Hamiltonian 
reads, 

a a 

+ JIP J2 Sa • Saa + jf P Yl Ma ' S » a 

a a 

+ + jf r T M a ] ■ (s lr + S rl ), (12) 

a=l,r 

where J?(D ) = #(£>„) = jT (jDq) = 

f^f , and J^(D ) = yll-jffipfiPo), and P = 

S a =i r (^ T " T ° +^ SaSa ). Here again, the Kondo physics 
falls under the category 25 = n. It is then sufficient to 
write (and solve) second order poor-man scaling equa- 
tions. Having done it, we find the stable fixed points 
Jf , jf T — ► oo (i = 1, 2). These define the corresponding 
Kondo temperature 

T ^ cxp ("(V3 + DK + jr))- 

For N — 4, similar analysis shows that the TQD again 
manifests the fully screened Kondo resonance with the 
(oo,oo) fixed point. 

Summary: The novel achievements are: i) In compos- 
ite quantum dots, such as the TQD displayed in Fig. 
n the two-channel (left-right leads) Kondo Hamiltonian 
emerges in which the impurity is a real spin and the 
current is due solely to co-tunneling. The corresponding 
exchange constant Ji r is a relevant parameter: by tak- 
ing even and odd combinations, the system is mapped 



on an anisotropic two-channel Kondo problem where Ji r 
determines the degree of anisotropy. ii) RG equations 
© for the exchange constants are solved iJSJ and yield 
the flow diagram displayed in Fig|2 iii) Although the 
generic two-channel Kondo fixed-point is not achievable 
in the strong coupling limit, inspecting the conductance 
G(V ga ,T) as function of temperature (Fig|3$ and gate 
voltage (Fig^J suggests an experimentally controllable 
detection of its precursor in the weak and intermediate 
coupling regimes, iv) Analysis of the Kondo effect in 
cases of higher spin degeneracy of the dot ground state is 
carried out in relation with dynamical symmetries. The 
Kondo physics remains that of a fully screened impurity 
n = 25, and the corresponding Kondo temperatures are 
calculated, v) Electron occupation is intimately related 
to the nature of the Kondo physics (n > 25 for N = 3 
and n — 25 for N = 2,4). vi) One remark regarding 
non-linear response is in order. If a quantum dot in the 
two channel Kondo regime is subject to a finite DC bias 
there is a peculiar effect of oscillatory current response 
[l2^ . The TQD in series discussed here might then be 
used to test it experimentally 
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